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Abstract 

The operational approach to the measurement of phase studied by Noh, 
Fougeres and Mandel is applied to the measurement of the state of polar- 
ization of fully polarized light. Operational counterparts of the quantum 
Stokes parameters are introduced and their fluctuations are examined. It is 
shown that, if the polarized field is weak, the measured fluctuations are in- 
fluenced not only by the quantum properties of the source field but also that 
of the measurement. This character is reflected on the measured probability 
distributions of the parameters of polarization, which are also investigated 
independently for the fully polarized coherent states and the Fock states as 
the initial field strength is varied. Finally, connection between the operational 
approach to the measurement of polarization and the su{2) interferometry is 
examined. 

PACS numbers: 42.50, 42.50.B, 42.50.D, 85.60. G 
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I. INTRODUCTION 



The idea of the operational approach as an experimental technique based on photon count- 
ing in the measurement of quantum phase fluctuations has been suggested in the late 80s 
by Barnett and PeggB in the context of a measured phase operator using certain homo- 
dyne experiments and more recently was formulated in detail by Noh, Fougeres and Mandel 
(NFM)i. The operational phase measurement is based on using N-port quantum homodyne 
detectors of which the analogy with classical homodyne approach is based on the purpose of 
extracting information about the phase between two initial fields by performing a complete 
set of photocount measurements between the components of the field. This procedure of 
obtaining the phase information between two fields depends on the particular experimen- 
tal scheme through its classical analogy of relating the relative photocount measurements 
to certain cosine and sine functions of the relative phase. Since, through this suggested 
analogy, different quantum measurement schemes would correspond to different classical 
ones, the information extracted for the relative phase is expected to be different for dif- 
ferent experimental schemes. Indeed this point has been demonstrated in the formulation 
of the operational phase measurement by NFM by starting with two different classical and 
quantum measurement schemes where one measurement used two-port homodyne detection 
whereas the second one used four-port homodyne detection!!!. The two port measurement 
yields either the cosine or the sine information about the phase failing to give the full phase 
information. In the four port scheme the simultaneous measurements were made possible 
by well-defined trigonometric operators of the relative phase where the full information on 
the phase and its fluctuations can be extracted. On the other hand a comparison of NFM's 
operational approach with the operational approach introduced by Vogel and SchleicbH has 
been compared by LynchU who found agreement between the two operational schemes. 
Another particularly important part of this scheme dependence manifests itself in the weak- 
field measurements in which the quantized nature of light as well as of the detectors becomes 
crucial when the homodyne detectors have a relatively high probability of registering a few 
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or null photocounts within the measurement time interval T. This being the case for a 
single set of measurements, one considers an ensemble of repeated measurements under the 
same initial conditions. There, each repeated measurement would have generally different 
but equally acceptable configurations of detected photons and one has to make a distinction 
between the outcome of a single measurement from the average outcome of a collection of 
such repeated measurements under otherwise the same conditions. Despite the fact that 
the experimental verification of the operational NFM approach was successfully made by 
the same groupi, the appearance of the discrete outcomes in the phase measurements in 
their scheme was subject to long and heavy discussions@~0. In this work we suggest another 
application of their approach to the operational measurement of the state of polarization of 
a fully polarized source. In an earlier publication,^ we investigated a particular extension 
of NFM's operational approach to the measurement of the Stokes parameters of a fully 
polarized weak coherent light. In this work, we will extend this formalism introduced in 
[11] to a more general framework by including the calculations for the measured probability 
distributions of the polarization fluctuations and also examine the case throughly when the 
initial field is a fully polarized Fock state. 

We start with a brief outline of the operational approach to the measurement of polarization 
fluctuations when the polarized field is given in a classical as well as a quantum state. In 
Sec. II we present the general formalism of calculating the polarization fluctuations and their 
corresponding probability distributions. Sections II. A and II. B are devoted to the specific 
calculations corresponding to two different fully polarized initial quantum states of the field 
as coherent and Fock states respectively. Section II. C is devoted to the connections between 
the operational approach and the su(2) interferometry. 

Classically, the state of polarization of a fully polarized monochromatic field, Ei = 
6i cos {ojt + 5i), where i — 1,2 are the polarization indices of two preselected orthog- 
onal polarization eigenmodes, can be manifestly described by four Stokes parameters 
s m (m = 0,1,2,3) asEHl 
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*o= \{(El) + {(El)) , 



*i= \{(El) ~ ((ED) , 

(1) 

s 2 = {(El) (ED) 112 cos0, 

, 3 = ((^ 2 )(i? 2 2 )) 1/2 sm^ 

where = 62 — 5\ is the optical (temporal) phase and U = (Ef) is the intensity of the 
corresponding ith component ( i = 1,2). We now describe an experimental setup based on 
a set of photocount measurements for the purpose of investigating the fluctuations in the 
measurement of the classical Stokes parameters in Eqs (fj) and their corresponding quantum 
counterparts. 

A. The Classical Measurement Scheme 

Within the operational approach, it is possible to measure all classical Stokes parameters 
in terms of the various components of the intensity. The experimental scheme is shown in 
Fig. (1). The initial field enters the setup through the %50 — %50 beam splitter BS±. One 
of the output beams of BSi is sent to a polarizing beam splitter PBS\ which defines a 
reference frame 1, 2 for the relative angular orientation of all other polarizing beam splitters. 
The other arm of the beam leaving BS\ is sent to BS2 as an input, leading to the second 
part of the experiment where the simultaneous measurements of cos and sin are realized 
independently from the first part. PBS2 is aligned at a 45° angle with respect to the 
reference frame selected by PBS\. The intensities measured at the detectors D 3 and 
yield the measured values of cos0 and its moments. For the sirup measurement, the phase 
of the remaining arm of the field is shifted by 7r/2 via a quarter wave-plate A/4. The field 
is then sent to PBS3 which is aligned in parallel to PBS2- A simple calculation shows that 
the classical intensities measured at all detectors Dj ,(i = 1, .., 6) are given by 
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h = \{E\) , 
l i = \{Ef) , 

h = \[(El) + (El) + 2 v /(£ 1 2 >(£ 2 2 >cos0] , 

(2) 

h = \[{El) + {ED - 2^(^)^1)008 0] , 

h = \[{El) + {ED + 2yl W) {ED sm0] , 

h = \[{El) + {ED ~ 2^<£ 2 )<£ 2 2 )sin0] . 
Eqs ([!]) and (Q) imply that the classical Stokes parameters can be extracted operationally 
by measuring all field intensities lj (i = 1, ..,6). In terms of these intensities, the Stokes 
parameters are simply given by 

s = (h + h) , si = (h - h) , 

(3) 

S2 = (h - h) , S3 = (h - h) ■ 

In [11], we parameterized the polarized field in terms of the functions 

cosO = Si/s , sin6> = \J Sq — s 2 /s , 

(4) 

cos0 = s 2 /\Jsl + si , sin0 = s 3 /-y/s2 + s 2 . 

This particular choice of the parameters above proves to be very convenient in the quan- 
tum operational measurements on fully polarized light. They also lead us naturally to the 
Poincare's geometric interpretation of polarization.&^l Here 6 and are physical parame- 
ters as shown in Fig. (2). Their values are directly connected with the ellipsometry of the 
polarized field. From now on we will adopt this parameterization and study the polarization 
fluctuations in terms of the fluctuations of these observables. 
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B. The Quantum Measurement Scheme 



The classical description above is adequate when the field intensity is sufficiently high. The 
vacuum fields, which are not present in the classical approach, are necessary for the correct 
quantum description of the apparatus as well as the field observables. 

In Fig. (1), the field operators di, d 2 at the output of PBSi are related to the input field 
components 01,02 as0 

di = (rai +tv{ 1) ), d 2 = -y= (ra 2 + tv { 2 ) ) (5) 

where r = i/y/2 and t = l/y2 are the field reflection and transmission coefficients, and 
(j = 1,2) are the polarized vacuum fields entering through the vacuum port of BSi. 
If the measurement scheme in Fig. (1) is extended to include the photodetectors Di (i = 
3, 4, 5, 6), then measurements can be made compatible with a proper quantum treatment 
of all fields. The output fields of PBS 2 and PBS 3 at D 3 , D4, D 5 , D 6 are given by 

d 3 = ^[(iro! + r 2 v{ 1] + tv^) + {tra 2 + r 2 v { 2 1] + tv { 2 ] )\ 



d 4 = 4j[-(*roi + r 2 ^ + tvf ] ) + (trd 2 + r 2 v { 2 ] + tv { 2 ] )} 



d 5 = -^W 2 di + trv[ 1] + rv{ 2) ) + (t 2 d 2 + trv { 2 1] + rv { 2 ] )} 



(6) 



4 = 1 [-i{t 2 ai + trv { i ] + rv {2) ) + {fd 2 + trv { 2 ] + rv { 2 ] )} . 

In connection with their classical counterparts in Eqs(|3|), we are now at a point to suggest 
the quantum Stokes parameters for the field operators dj within this operational approach 
in terms of the observable photon number operators rij = d\ di as 

S = n x + n 2 , S x = hi - h 2 , 

(7) 

S 2 = ^3 - n A , S 3 = h 5 - h 6 . 
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In Eqs(6), all field operators commute as a manifestation of the vacuum fields. Hence in 
Eqs (|7|) we have [Ej, Ej] =0 (i ^ j) and all photon number operators can be simultaneously 
measured at the detectors Di (i — 1, . . . , 6). As a result, Eqs (0) are compatible with their 
classical counterparts in Eqs (|3|). This property of the Ej , (i — . . . 3) operators allows us 
to further suggest an extension Eqs (f|) to their operator counterparts as 

C e = E^Ei , S e = (1-(% 2 ) 1/2 and 

(8) 

(% = E 2 (EI + Sl)-V2 , ^ = S 3 (El + El)-V2 . 

Cq, Sq as well as S$ are well-defined and compatible quantum observables. They commute 
with each other and satisfy the operator relations Cg + Sg = 1 , C*| + S| = 1 and, as the 
result, they can be measured simultaneously. 

One of the benefits of adopting Eqs ([?[) and (|8|) is that all measurements are now based on 
pure photon counting depending on the measured photocounts at the detectors Di , {i = 
1, . . . 6), and hence they do not involve any temporal interference effects. This is an advan- 
tage of the operational measurement, which will be transparent later in our discussion of 
the weak-field limit. 

Eqs®, hereinafter referred to as the operational quantum Stokes parameters (OQSP), are 
the most convenient choice for Ej (i = 0, ..,3) befitting the purpose of the photocount 
measurement scheme of Fig. ([]]). All operators in (||) are now compatible with the classical 
variables of Eqs (Q) as long as the measurements of the E 2 and E3 operators do not yield 
zero simultaneously. 

II. THE MEASUREMENT OF POLARIZATION FLUCTUATIONS IN WEAK 

FIELDS 

The operational approach as applied to the polarization measurement of a fully polarized 
and weak initial field is based on individual detections of single photons where the quantum 
nature of the field as well as of the detection mechanisms is dominant. The influence of the 
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direct quantum homodyne detection on the statistics of a quantum measurement has been 
examined by MandelEl, Kelley and KleinerO as well as by Glaubei§ and expressed for in 
the form of a combined quantum probability distribution 

N 

V({n j }) = l[ : (S^P expi-^/njl : . (9) 

where : : accounts for the normal ordering of the operators d\, di inside and d\ di corresponds 
to the photon number operator. Throughout the calculations the measurement time interval 
will be assumed to be much smaller than the coherence time (which is naturally satisfied for 
a monochromatic field) and much larger than the inverse of the oscillation frequency of the 
field. Under these conditions it is possible to consider the simplest case when the photocount 
measurement at the detectors are time translationally invariant and linearly dependent on 
the measurement time interval T. 

Including the quantum effects of the homodyne detection in Eq. (^), an individual measure- 
ment of an arbitrary field operator f{{fij}) yields the measured value 

(/({%}))= AT £ f{{n 3 })Tr{pV{{n 3 })} j = l,2,...M (10) 

{%*} 

where the trace is considered over the complete set of states in the density matrix of the 
initial field p= \ip) in in (ip\. 

With Eqs ([D and flT0|) representing a general scheme of measurement in the operational 
approach, we now consider for /({%}) the operators of {fij} , (j — 1, 2 or 3, 4, 5, 6) 

£ e (x) = (C e + i S e ) x , and £^(x) = (Q + i S^) x where i£R. (11) 

In the construction of Cg,Se and C^.S^ pairs in Eqs(H), the compatibility conditions 
= of the OQSP ensure that = 1 and ||£^(x)|| = 1, hence Sg(x) and 

S^(x) are unitary operators for all x £ R. According to the procedure outlined in the 
context of Eqs (P) and fllOD , the measurements of these operators yield 

E + (12) 



and 

(4M> =K E [ ^-"*> + i <"°-"°> ]■ (P ({ .,, }) > (13) 

yj (n 3 - n A y + (n 5 - n 6 )^ 

where (P({n,») = Tr{\^) in i n (i>\V({nj})}. In Eq. (0) {n,} = (m,n 2 ) and in Eq. ([TJ 
{ n j} = ( n 3' n 4' n 5' n 6)- Clearly, Eq. ([12]) is well defined if rii,n 2 are not simultaneously zero 
and, similarly, Eq. flilf ) is well defined if (n 3 — n 4 ) and (n 5 — n 6 ) are not simultaneously 
zero in the respective summations above. The idea of the elimination of the configurations 
n i = n 2 = and (n 3 = n 4 ), (n 5 = n 6 ) from the statistical weight has been introduced as 
a crucial element of the operational approach!!!!'! in the implementation of the statistical 
averages. The effective weight of such configurations becomes non-negligible particularly in 
the case when the initial field strength is sufficiently weak where the probability of receiving 
zero number of photons within the detector's measurement time interval T is finite. For 
instance, the weight of observing zero photons simultaneously at the detectors D\^D 2 is 
given by (V(0, 0)). The result of such a null measurement is inconclusive in the calculation of 
the averages in Eq. flT2|). Similarly, n 3 = n 4 and n 5 = n 6 yield additional inconclusive results 
in the measurement on (S^(x)) in Eq. (|T3D. The measured averages are then normalized by 
excluding the total statistical weight of these inconclusive configurations from the integrated 
probability. For strong fields, the weight of such ambiguous outcomes is smaller and in the 
classical field limit there is no contribution from such terms, viz., Me = M^ = 1. In the 
measurement of the temporal phase the individual fluctuations of these weak components as 
well as the fluctuations in the relative number of photons can be strong due to the absence of 
a classical reference source (i.e. a strong local oscillator). Hence the normalization technique 
introduced by NFM proves to be essential for any operational measurement based on phase 
and thus also for our approach here. 

More explicitly, this normalization procedure amounts toi 

Mi 1 = 1-(V(0,0)) (14) 

and 



A^, 1 = 1 - (^( n 3 = n 4 = n,n 5 = n 6 = m)) 



(15) 



in Eqs flT2]) and flTB|). The observed unitarity conditions of £q{x) and £^{x) suggest that one 
can associate a classical random variable e lxe and e lx< ^ respecting the probability distributions 
P(9) and P((f>) such that! 



(S (x))= / d6e lxe P(6) , and (^(x)> 
Jo 



.e ia ^P(</>) . 



(16) 



The probability distributions can then be obtained by the inverse Fourier transformations 
of © by 

/CO Ann , A A - /"CO /-7/t> a 

_ e -*' {<£(*)) + (**(-*))} , and P(0) = J — e -*+(E+(x)) 



(17) 



with d6P(e) = JZ ir d<l>P^) = l. 

Defining two auxiliary functions of {n}j by 

_ 1 /2 A /nin 2 



tan 



and 



W 



tan 



1/ ^5 ~ ^6 



(18) 

n x - n 2 y v n 3 - 

where {n} = (m, n 2 ) and {n} = (n 3 , n 4 , n 5 , n 6 ) for 6 1 and respectively, and using Eqs (|l2 
and (|13|), the moments for a generalized initial state \ip)i n read 

(4(x)) = J\f E' 2 e**°w (V({ nj })} where 



"2 



(ni— n2)+«2 v /ni ri2 
ni+n 2 



(19) 



x : (a[ai) ni (a5,a 2 )™ 2 exp[-l/2(a[a 1 + a^a 2 )] : 



2 "i+"2 ni !n 2 ! 

t. 



and 



<4(a;; <?„)>= E' r 



n3, n.4, n 5 ,n 6 



e w*w <J )> where 



n3,n4,n 5 ,n6 



(n3-n 4 )+i(n 5 -n 6 ) 
^ («3-n4) 2 +("5-"6) 2 



1 

8 n 3 +n 4 + n 5 +n 6 „ 3 l„ 5 ! ng I 



X 



: (a\ + 4)™ 3 (a, + a 2 )™ 3 (-a{ + a 2 )™ 4 (-a x + a 2 )™ 4 



x (-vial + 4) ns + a 2 )™ 5 (ia5 + a\) n6 {-ia x + a 2 ) m 
x exp[-l/2(a\a 1 + 4a 2 )] : \ip) in . 



(20) 
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with 5o = 62 — 8\ implicitly described in Eq. (pO| ) as the relative temporal phase between the 
components of the initial field. The primes on the summations in Eqs (|T9| ) and ([H]) now 
indicate that the summations are performed by excluding those configurations for which the 
outcome is inconclusive. 

All moments are now determined once the initial components {fix), (^2) and the relative 
temporal phase 5q of the inclusive fields 01, (12 are known. In our calculations the initial field 
parameters are chosen as the ratio of the photon numbers rj = {n\) / (^2), the total number 
of photons (So) = (ni) + (^2), and the relative temporal phase 5q. 

The credibility of the results obtained from the quantum operational approach crucially de- 
pends on the understanding of the influence of the quantum detectors on the final statistics. 
As pointed out in Refs. [2,3,6,9] that another essential element of the operational approach 
is the construction of an ensemble from a long series of such single operational measure- 
ments. The final physical results are then obtained by averaging the outcomes of single 
measurements over the created ensemble. Based on this prescription, we must now con- 
struct a physical ensemble of measured configurations in the calculations of the moments 
as well as the probability distributions in Eqs flT3) and fl2"0|). The response of the quantum 
detectors to the incoming photons in the creation of the photocurrent is a random process 
which obeys the Poisson statistics in Eq. (H).0 As the photoelectrons are emitted at random 
times respecting this statistics, the information regarding the initial temporal phase 5q of 
the incoming photons is modified and each repeated measurement is equivalent to super- 
posing a random phase shift A on 5q. Hence the process of repeated measurements creates 
an ensemble of temporal phase configurations 5q + A with A being uniformly distributed 
over the available range. Since we consider in our calculations that the measurement time 
interval T is considerably larger that the coherence time, the available range for A is the 
entire 2tt range. Hence the average over the created ensemble corresponds to an averaging 
over a uniform distribution of A. It is clear from Eq. ( |i9|) that the moments (Sg(x)) are 
independent from 5$; hence they will also be independent of A. This implies that a uni- 
form average over A does not influence the measured moments (Sg(x)) and the probability 
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distribution for P{6) is given by 

dx 



— e- ix0 {(S e (x)) + (£e(-x))\ , < Q < n 



(21) 



On the other hand, the moments (S^(x)) depend on the temporal phase Sq and before the A 
average, Sq dependence must be replaced by Sq + A. This produces, at each measurement, the 
conditioned moments (£<f>(x; S + A)) and, following Ref. [6], their conditional probability 
distribution P(0, Sq] A) is given by 



/OO A <J> 
^(^(x;5 + A))e—^- A ) 
-oo ZTT 



Therefore, the ensemble averaged probability distribution is 

-* dA 



-7T < < 7T 



(22) 



P(<>,6 ) = y ^ — P(0,5 O ;A) 



(23) 



After a short calculation using Eqs ( |T9D and ( |20| ) in Eqs. ( |21| - p3| ), the probability distributions 
P{9) and P(0, Sq) can be expressed by 



P{9)=N £ 8(0 -0 n ) (V({ nj })) 



(24) 



{nj} 



and 



p(0, 5 ) =^J2 CP(M, to - 0„ + 0)) • 

{rij} 



(25) 

where the last term in Eq. (|25|) is obtained by using Eqs ( p0|) in Eq. ([22]) . 
On the other hand, the detectors' influence on the measured statistics can only be understood 
if the measured moments and probability distributions are compared with those without the 
detectors' influence. For this purpose and, following Refs. [2,3], we define the theoretically 
inferred values of the 9 and moments as (Sg(x)) and (S^(x)) where 

-hi — n 2 + i 2 y/hi h 2 - 



(£${x)) = m M 



hi + h 2 



(26) 



and 



(h 3 -h 4 )+i (n 5 - rag) 
n 3 - n 4 ) 2 + (n 5 - ra 6 ) : 



(27) 
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where : : stands for the normal ordering of the field and vacuum operators inside. 
We calculate the probability distributions P(9), where < 9 < n, and P(<P,Sq), where 
—7i < <p < 7T numerically using Eqs (^) and (p5|). Since cos 9 is single valued in the 9 
range considered, we will only need to examine the fluctuations in the Cq operator. On 
the other hand, in the range considered both and operators will be necessary. In 
our calculations, the summations over infinite range of {n}/s are truncated at {n}™ ax = 20 
for all j which naturally restrict the accuracy of the results to sufficiently weak initial 
fields. The measured moments and the probability distributions are then compared with 
the theoretically inferred ones by using Eqs (^) and (|27| ) . 



A. Calculations for a fully polarized quantum coherent field: 



Let us now assume that the initial field is in a fully polarized quantum coherent state 
\i>)m — \ a i, a 2), with the parameters given by ctj = \aj\ e lSj where |<x,-| 2 and 5j (j = 1,2) 
are the average number of photons and the coherent temporal phase of the j'th component 
respectively. The relative temporal phase is given, as before, by 5$ = 62 — Si. From Eqs ([B5]) 
and (|20|) , the measured moments in this state are given by 

(S e (x)) = Me Y! nu 



,n 2 



(ni-ra 2 )+t2 v /?ii ri2 1 x 
ni+n 2 



x |ai| 2ni \a 2 \ 2n2 exp{-i(|ai| 2 + |a 2 | 2 )}/2 ni+n2 m!n 2 ! 



and 



(£t(x;5 ))= A/; Ei 



(n 3 -n,4)+i(n 5 -n 6 ) 



■y/ {nz~ni) 2 +(n 5 -ne,) 2 



(28) 



X Oil + «2 



|2n 3 



«1 + a 2 



2n4 



%Qi\ + «2 \i0t\ + «2 



(29) 



xexp{-i(|«i| 2 + |a 2 | 2 )}/8 n3+n4+n5+n6 n 3 !n 4 !n 5 !n 6 ! , 
where {n} = (n 3 , m, n^, uq). For the specific initial polarized coherent state considered, 
using Eqs(|l4D and ([Top , the normalizations are given by 
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Mi 1 = 1 - exp{-I(| ai | a + |a 2 | 2 )} (30) 

and, denning (3 = 2\ai\ |a 2 |/(|«i| 2 + |ck2 | 2 ) where (3 < 1, 

K r-x _ i _ v ^M!±M^A 2{n+m) exp{-\{\ ai \+\ a2 \ 2 )} 
JV <t> ~ 1 ^n,m{ 8 ) (n!) 2 (m!) 2 ' 

l-f3 2 cos 2 5 ]" [l -(3 2 sin 2 5 ] m 
We will first examine the P(9) distribution. Using Eq. ( |24] ) the calculation of P(9) yields 

P(9)=M e £ 6(9-9 {n} ) ' 1 ' ^ , ex P {--(K| 2 + |a 2 | 2 )} (32) 

where 9{ n } is defined by the first expression in Eqs (|T8D . For sufficiently weak fields, i.e., 
(S ) ^ I? each detector measures null or a very few number of photons. This implies that 
in Eqs (|32]) it is sufficient to restrict the summation over {nj} , (j = 1,2) to a few terms. 
For instance, let us consider {nj} = 0, 1. Then including only the first-order terms in the 
average total photon number, Eq. ( p2[) can be approximately expressed in the weak-field 
limit by P w (9) in the form 

P w {9) = Af w (S ) { — 5(9) + 5(9 - tt) ) (33) 
[ 1 + T)~ x 1 + 1] J 

where (So) = (|ai| 2 + |a 2 | 2 )/2 is the total average photon number deduced from the mea- 
surements at the detectors Dx,D 2 , r\ = |ai| 2 /|a 2 | 2 and M w = (So) -1 so that Jq d9 P(9) = 1. 
From Eq. (|33| ) we find that 

(cos9) w = r d9 cos9P(9) = , (cos 2 9) w = l (34) 

Jo r/ + l 



Clearly, (cos9) w in Eq. (|34]) is consistent with the theoretically inferred values calculated 
from Eq. (|26|), [i.e. (cosfl)™ = (Cg)]- In the initial polarized coherent state the theoretically 
inferred moments are given by 

isUx)) = vn- l + i27 r' 2 y = e -tan- 2 ^/(.-D (35) 

L 77 + 1 J 

which respect a nonfluctuating distribution. Equation (^) is also consistent with the clas- 
sical calculations using Eq. @. However for the second moments we obtain 
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{{Clf) = (^-|) 2 ^ (cos 2 9) w . (36) 

The 6 1 distribution in Eq. (^) is plotted in Fig. 3 below for r\ = 1.0,0.5 and (So) = 
0.1, 1.0,5.0, 10.0. The first observation in Fig. 3(a) is that at i] = 1.0 the probability distri- 
bution is symmetrically centered around 9 = ir/2. In the weak-field limit P(9) is peaked at 
9 = 0, 7r. As the field strength is sufficiently increased, the central peak at 9 = tt/2 gradually 
develops as all other peaks are suppressed. The average of the cos# within the full range 
< 9 < 7i is zero as it would also be expected from the theoretically inferred moments 
in Eq. (|35|) . For rj ^ 1, the measured P{9) is plotted in Fig.3(b). The delta functions in 
Eq. ( p2|) are numerically simulated by sharp Lorentzians, hence they acquire a finite width in 
Figs. 3(a) and 3(b). On the other hand, using the equation in (f21|), the inferred probability 
distribution P I {6) can be found as P I (9) = did - cos" 1 ^ - l)/(r} + l) 



A similar calculation can also be done for the P(4>) distribution by making use of Eqs ([2^). 
(f22|) and (|23|). After some calculation using the normalization procedure leading to Eq. ( p3|) 
we find that 

r(,<p,0 ) - l^ {n . } ,(,-=3,4,5,6) n 3 \n 4 \n 5 \n 6 \ { 8 J 

X[l + P COS(5 " <P - <P{n})Y* [1-/3 COS(5 ~ <t> ~ 0{n})]"' 4 ( 37 ) 

x[l + /3 sm(5 - - M )] n5 [1 -13 sin(5 - - {n} )]" 6 

which is, not surprisingly, the same distribution obtained by NFM in Ref. [6] in a slightly 
different context. The weak-field limit of Eq. (|37|) has also been studied in Ref. [6] to where 
we refer the reader for additional details. The numerically calculated Eq. (B71) is plotted in 
Fig. 4 for rj = 1.0, 0.5, So = 0.1, 1.0, 5.0, 10.0 and for Sq = 0. The first observation we make 
here is that P(<f>, Sq) is almost independent from 7] but strongly dependent on the strength of 
the initial field. As the field strength increases, the fluctuations decrease and the distribution 
becomes gradually narrower. On the other hand, using Eq. (|27|), the theoretically inferred 
moments are calculated as 

= e**" . (38) 
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Hence, the theoretically inferred distribution is also non-fluctuating given by P(4>, 8q) = 
8 (<f> — 8q). The operational averages for <ft as well as the probability distributions are strongly 
peaked in the strong-field limit and they have the tendency to approach to the theoretically 
inferred moments and the delta function-like probability distributions respectively. On the 
other hand, the operational approach predicts large deviations of the measurement from 
the theoretical values in the weak-field limit. In order to understand the influence of the 
photodetection particularly in the weak- field limit, we now examine the second order fluc- 
tuations in the measured moments of the 9 and related operators. 

1. The measured fluctuations in polarization 

Once the moments in Eqs and (^) are defined, the measured moments of the cosine 
and sine operators of 9 and can be found. The same moments can also be found by 
the use of the probability distributions in Eqs (|32l) and fl37|) . Here the weak-field limit is 
particularly interesting and it can also be examined analytically. We start our analysis of 
the fluctuations by reminding that, since < 9 < tt, we will be confined to the measured 
fluctuations in the Cg operator. In the weak-field limit (keeping only the leading term in 
the total field strength) 

<<%> * ^ , (<% 2 > * 1 , (39) 
we find the dispersion D{9) as 

D{9) = ^{{Ce-{Ce)Y)^-^—<l (40) 

rj + 1 

The dependence of D(9) for 77 = 1.0, 0.5, 0.1 on the total field strength is shown in Fig. 5. We 
now shift our attention to the measured fluctuations in the <fr distribution. Since — tt < <p < tt 
we need to consider here both the cosine and the sine moments. Considering first the weak- 
field limit, and keeping only linear terms in the total field strength, have, 

liM^, (S^'^^ (Cj> = (Sj) = i (41) 
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In this case we define the dispersion D(<f)) as 



m = V((^-(^» 2 > + ((^-(^» 2 > = ^ - > ^ ( 42 ) 

The dependence of D(<f>) on the total field strength is plotted for rj = 1.0, 0.5, 0.1 in Fig. 6. 
On the other hand, for both the 9 and the related moments the theoretically inferred 
fluctuations vanish in the coherent state (i.e. D{9) = -D(0) 7 = 0). The measured fluctu- 
ations differ significantly from the theoretically inferred ones as the strength of the initial 
field becomes weak. These deviations in Figs. 5 and 6 from the theoretically inferred values 
arise from the nature of the photodetection of weak fields and the normalization of the prob- 
ability weight after discarding the inconclusive data. We find that the results for the fully 
polarized coherent field and the differences between the measured and inferred fluctuations 
in the weak-field limit closely relate to the results obtained by NFM.i'i 
The 7] dependence of the fluctuations in Figs. 5 and 6 implies that r\ can be used as a 
parameter in the measurement to search for an optimum orientation of the set-up in Fig. 1 by 
rotating the reference axes 1, 2 around the initial field direction. Note that this corresponds 
to a solid rotation of the entire setup since the relative orientation of each polarizing beam 
splitter with respect to PBS\ is fixed. By this operation the angle between the polarization 
axes of the setup and the mean principle axes of the polarization ellipse of the initial field can 
be changed. Let us suppose that a w and a 20 are initial coherent state parameters defined 
with respect to some fixed orientation 1o,2q of PBSi and given by ajo = \ctjo \ e t5j0 . If the 
principle axes 1, 2 are rotated by an angle 7 with respect to lo, 2o, the initial coherent field 
parameters a,\,a,2 are effectively rotated by the same angle with respect to a 10 and a 2 o- in 
particular, the average number of photons (77,1) and (n 2 ) measured at D\ and D 2 are given 
by 

1 2 1 2 

(ni) = - ctio cos 7 - a 2 o sin 7 , and (n 2 ) = - a 2 o cos 7 + a w sin 7 (43) 

Since 7 is arbitrary, we can use it to tune r\ = (ni) / (n 2 ) in order to find whether an optimum 
orientation of the setup exists such that both 9 and related measurements (or whatever 
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other observables are examined) can be improved simultaneously. The measured rj at the 

detectors D 12 is then a function of 7 and is represented in terms of the initially fixed 

Vo — l a io| 2 /|«2o| 2 and the relative phase (£20 ~ £lo) as 

/ s _ (VVo ~ tan7) 2 + ^rfc tan 7 sin 2 (S 20 - 5 w )/2 
(y/Vo + tan7) 2 - v /r^ tan 7 sin 2 (5 20 - 5 10 )/2 

The second order fluctuations represented by D(9) and D{(j>) are plotted in the weak- field 

limit as a function of 7 for < 7 < tt/2 and for (So) = 1 in Fig.7 and for (So) = 9 in 

Fig. 8. The figures imply that such an optimum orientation to simultaneously minimize the 

fluctuations D{6) and D(<f>) for a fixed value of 7, 5 2 o — ^10 and (S ) does not exist. Hence, 

depending on the measured observable, one has to engineer such optimum configurations for 

each measurement independently. 



B. Calculations for a fully polarized Fock state: 

Now let us assume that the initial field is given by the fully polarized photon number state 
|V>)i„ = \M)^ as 

m /m\ 1/2 

\M)* = Mm E e*° (M " m) \m , M - m) (45) 

m=0 \ m J 

where 0o is a temporal phase between the polarization components, Mm — 2~ M / 2 , is the 
binomial coefficient and \m, M—m) = \m)<g>\M—m) describes the relative number of photons 
in each component in reference to a particular choice of predefined axes of polarization. It 
is possible to see that the field operators for the Fock state 

= -j= (o! + e"** a 2 ) , and a\ o = -±= (a{ + e i<t>0 a\) (46) 

satisfy 

IM)^ = VM\M- I)*, , and a\ \M)$ = VM+l\M + 1)^ (47) 

The temporal phase factor 0o determines the ellipticity of the polarization. If 0o — 0, tt, 
Eq. ( |45"D provides the basis for linear polarization. For <p = ±7r/2 left and right circularly 
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polarized states are obtained. For arbitrary 4> left and right elliptically polarized Fock states 
can be produced. Restricting <f) within the range < 4>q < ir, the left and right elliptically 
polarized states are realized respectively by |M}^ and \M)^ _ n ) with the respective field 
operators d^ ,a^ o ; d^-^), oL _ % y The second pair of field operators are then found by 
making the change 0o ~~ * <Po ~ ft m Eqs (pUf). The field angular momentum operator is given 
by C z = (dj, — a,^_ n d^-^) and C z is diagonal in \M)<f, with eigenvalue M. 
Using Eqs ((46|) and (0), Eq. ( fj5| ) can be written as 



|M> *» = (^?^ m i°-°> = 7M < a *.' M i°> < 48 > 



where |0) is the vacuum state for a^ as well as for d\ and, d 2 ; (i.e., |0) = 1 , 0) ) . In what 
follows, the full range —it < <p < it will be considered. In fact, Eq (pE5|) is an example in a 
class of fully polarized Fock states corresponding to rj = (h\) / (fi2) = 1 where (ni) and (77-2) 
describe average number photons in individual polarization modes. For Eq. (^) we have 
(fii) = (fiz) = M/2. If a rotation parameter 7 is introduced [for instance, like in Eq. ( |43|) 
for the coherent state] in the field space by 

d^ 0j7 = (cos 7 ax + e~ l ^° sin7a 2 ) , and a^, = (cos7a{ + e 1 ^ sin 702) , (49) 

in terms of the new field operators a<^ 0i7 , a^ o 7 the field operators of the initial Fock state 
in Eq. ([IS]) would be obtained when 7 = 7r/4 in Eq. (^5|). This implies that the Fock state 
|M)0 O)7 created by Eq. ([4^) is realized effectively by a (7 — 7r/4) degree rotation of the Fock 
state in Eq. (f45|) with |M)^ 0i7 being 

1 M ( M\ 1/2 

|M) 0Oi7 = -= [d\^) M |0) = £ (cos 7 r(e^sin7) (M - m) KM-m). (50) 

VJM! m=0 \ m / 



Eq. ( |50"D for a fixed 0o now describes a fully polarized generalized Fock state with an arbitrary 
ratio of photon numbers 77(7) = cot 2 7 between the polarization components. 
In comparison with the coherent initial field, a considerably more tedious work is involved in 
the numerical calculations of both measured moments. In the general fully polarized Fock 
state given by Eq. fl50|) , Eqs ([12]) and ([TJ) become 
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r(m - n 2 ) + i2y/n{ n 2 



ni,n 2 



and 



(^(x,0 o ))=A^E / M 



ni + n 2 



(n3~n4)+i(n 5 -ne) 



y/ (n 3 -ri4) 2 +(n 5 -na) 2 



2 n i+ n 2 n \ n \ y^—f 2 r+ Pr\p\ 



8™3+"4+«5+™6 n 3 ! n 4 ! n 5 \ n 6 ! 



9H-p r !„! } ( 51 ) 



(52) 



2 ai\ n 4 / .a[a 2 -4ii Nn: ' 
2^2/ V ajai+aj,a2 - 



• a|a2— alai \ n 6 



a[ai+a 2 <J2/ 



where in Eq. (|5T|) 



(53) 



Z^m=0 



M! 



COS 



7) 2m (sin7) 2(M_m) 



(m—n\—r)\ (M —m—ri2—p)\ 

with ni + n 2 < M in Eq. (^l|) and n < M where n = + + n 5 + n 6 in Eq. fl52"|) . In 
Eq. (|52"D, (£^(x,<f>o)) is understood in the same sense as (S^x, 8 Q )) in Eq. fl20"l) . 
The normalizations are determined as before by satisfying the condition (£g(0)) = 
(4(O,0 O )) = 1. 

The simplest analytic results can be obtained for the case M — 1 with 7 and 0o being free 
parameters. This corresponds for the initial state to 



|l)^ 0i7 = cos 7 |l,0) + e^ sin 7 |0, 1) 



(54) 



which is a fully polarized version of the split photon state in [2,6]. Using Eqs ( ^l| - |53|) we find 
for the moments 

Ma 



(^(x)) = ^[cos 2 7 + (-rsin 2 7 ] 



(55) 



where M e 1 = 1/2 and 

(£^x),<Po) = ^ { [1 + (-)* + i x + (-i)l + [1 - HI sm 2 7 cos 0o + \i x - (-i) x ) sin2 7 sin, 

(56) 
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where M^ 1 = 1/2. For the probability distributions P(6) and P((f), 0o) we use Eqs. (|2l] - [2~5"l) 
in the same spirit as we applied to the coherent initial state in Sec. II A. A simple calculation 
yields that 

P(6) = cos 2 7 5(0) + sin 2 7 (5(#-7r) (57) 

and 

P(0, O ) = 7T i 1 + sin2 ^ cos (^ " 40} ( 58 ) 

where the probability distributions are positive definite and properly normalized, i.e. 
Jq dO P(0) = dcj)P((j), 0o) = 1- At this point, a crucial limiting case in Eq. fl5"T|) and 
(|58|) needs to be mentioned. For 7 = 7r/4, Eqs (0) and (|58| ) describe the probability distri- 
butions of a fully polarized symmetric Fock state. For 7 = and tt/2 we have single mode 
photon Fock states |1,0) and |0, 1). The measured probability distributions in these states 
are P{6) = 5(6), P(0, O ) = 1/2tt for 7 = and, P(6) = 5(6 - tt), P(0, o ) = 1/2tt for 
7 = 7r/2; which, correctly describe the statistics of the single mode Fock state consistently 
with the theoretical expectations of a uniform distribution for P(4>, 4>o)- For all other 7, 
Eqs fl57|) and (|58D correctly describe the theoretical distributions for a general |M)^ 07 . This 
behaviour of the probability distributions can also be observed in Eq. fl3T| ) in the limits 1 <§C 77 
and f|< 1. The analytic calculations become exponentially harder for 2 < M. Nevertheless, 
explicit forms of the P(6) and P(0, 0o) can be given for a general M as 



m = A/i £ - * w ) {£ « } (59) 

where Eq. is used, and 



■ /v </> Z^{n} 8™3+™4+"5+"6 n3!n4!n 5 !n 6 ! ^r=0v 2 / r! £->p=0 \pj 

^3=0 ^ 3 J ^4=0 \e 4 ) ^5=0 ^^6=0 v^6 / 



^fc 3 =0 Vfc 3y ) Z^fc 4 =o ^fcj ^5=0 Vfc 5 y ^^6=0 Vfc 6 y ^ L J 

(cos^) 2m ^ e+k (sin 7) 5 

(60) 



x E™-n e-^-^w)^-*) _ -M _ (cos 7 ) 2m -^ ( s i nT )2(M-m)+l-fe 
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with and 0{ n } as given by Eqs. fll8D - The numerical calculations of Eqs. ( |59"D and ( |60D for 
linear polarization (e.g. O = 0), for i] = 1.0,0.5 [i.e., corresponding to 7 = 7r/4, tan _1 (\/2)] 
are presented in Figs. 9 and 10 for various values of M. Like in the coherent case, the 
temporal phase factor O in Eq ( |60| ) only shifts the distribution and does not play any role 
in the fluctuations. We now shift our attention to the second order fluctuations in the 9 and 
dependent moments. 



1. The measured fluctuations in polarization 

Similar to the coherent state example in Sec. II A, we can examine the 9 and dispersions 
in the weak-field limit in the range < 9 < 71 and — 7r < < 7r using the same observables 
as in Sec. II. A. 1. For M = 1 we have for 9 

((%>=cos2 7 , (C 2 e ) = l (61) 

hence 



Do = \[<Spe ~ (Ce)) 2 ) = sin2 7 = ^ (62) 



and for 



(^0) = ^ sin2 7 cos0o , 04,) = - sin 27 sin O , (Cj) = (Sj) = - (63) 



hence 



1 ™2 o„, _ ft V ^ 



d{4>) = V((^ - ^)Y) + m - (^)) 2 ) = y 1 - 4 sin 2 7 = yi - > ^- • (64) 

It is not an accident that the weak-field limit for the coherent state described in ( f£2p coincides 
with the Fock state calculations in ( p^) for M = 1 for all 77. The results of the numerical 
calculations of D{9) and D(<fi) as the initial number of photons is varied are shown for 
rj = 1.0, 0.1, 0.01 in Fig. 11 for D(9) and, Fig. 12 for D((p) corresponding to linear polarization 
(e.g. 0o = 0). Due to the large number of summations in Eq. ( |5"0"D calculations are considered 
within the range 1 < M < 10. 
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The fully polarized Fock state is a typical example where the correlations are present between 
the relative occupations of the polarization components. As pointed out in Ref. [2], this 
renders the physical interpretation of the theoretically inferred moments for the 0-related 
operators impossible. It appears that the operational approach here provides a scheme 
where the temporal phase distribution can be measured even if such correlations are present. 
We believe that the results obtained in Figs. 9 and 10 should be checked experimentally 
with the particular emphasis on the weak-field regimes, which we expect to provide further 
confirmation of NFM's operational scheme. The split-photon state discussed by NFM in 
[2,6] can also be interpreted as the weak- field limit of the polarized Fock state in Eq. ( f4~5|) 
corresponding to M = 1 where strong intensity correlations are present. For this state the 
inferred moments of the corresponding and are unphysical because of the fact that 
in Eq. (|27j) the denominator vanishes. To examine the theoretically inferred moments for 
a general M we use Eqs (|26"|) and ( p7|) in the Fock state ([45]) (we consider 7 = ir/A for 
simplicity) to calculate 

(M| : (ni -h 2 f : |M> 



<(<?/) ) 



(M| : (ni + n 2 ) x : \M) 



(65) 



where we find that 



{cD = {{cl) ) = 



(66) 



and 



(M\: (h 3 - n 4 ) + i (n 5 - n 6 ) : \M) (: (a\a 2 + a\a\) + i{a\a 2 — a\ai. 



) :> 



(M|: 



n 3 - n 4 ) 2 + (n 5 - n 6 ) 2 



:|M) 







d\ CL 2 (liCl 2 



(67) 



The vacuum fields do not contribute to the normal ordering and we also omitted the state 
label M in the second step of the expression. In order to calculate Eq. (|6Tf) we need (: 



\/fi\fi 2 :). In the presence of correlations (i.e. {n\fi 2 ) 7^ (h\)(n 2 )) we have 



{h\h 2 ) 
{hi){h 2 ) 



1 < 1 

M 



(68) 
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Hence, the correlation effects cannot be ignored if the initial Fock state contains a few 
number of photons. Furthermore, Eq. (|68| ) implies that, for M = 1 the denominator in 
Eq. ( |67D diverges at x = 2. It might therefore be suggested to consider that the comparison 
with the theoretically inferred moments with the measured operational ones is limited to the 
strong-field regime (1 <C M) where also consistency with the classical results are expected 
to hold. On the other hand, the denominator in Eq. (|67|) is not well-defined for values of x 
not equal to an even integer. Now let us assume for the moment that we are able to replace 
the denominator of Eq. (^7]) by ((dla^i^))^ 2 ■ One expects that if this replacement can be 
done, it can only be valid in the sufficiently strong-field limit where the correlations as well 
as fluctuations are expected to be negligible. With this replacement, Eq. (^7]) would yield 
cos 0o / qI\ sin 0o 2 , kqI\ 2 \ ■ 2 , 

(69) 

The inferred dispersion -D(0) 7 calculated from Eqs ( p5| ) is purely imaginary for all M which 
is an unphysical result. Hence, the replacement we made above, in order to make the de- 
nominator of Eq. ( pTI ) calculable, is unphysical for all M; thus, it cannot be done. Unlike the 
coherent state, the comparison with the theoretically inferred moments is made impossible 
by the presence of strong correlations. Therefore we are unable to examine the photode- 
tector effects in the weak-field limit in the operational measurement of the fully polarized 
Fock state using the standard formalism of theoretically inferred moments. The unphysical 
results we obtained for the inferred moments are not inherent to the quantum scheme. Even 
in the classical measurement scheme, there is no unique way of extracting the theoretically 
inferred moments when the relative phase or the relative intensity fluctuations are corre- 
lated. We refer the reader to Ref. [2] for a detailed discussion on this topic in the context 
of operational phase formalism. Nevertheless, we will suggest in the following subsection 
that for the Fock state in Eq. (|50|), or specifically for Eq. flilf), it is possible to find another 
measure to examine the photodetector effects in the weak-field limit by making use of the 
properties of the uncertainty relations. 
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C. The fully polarized Fock state and connections to the su(2) interferometry 



For symmetric distribution of photon numbers in the components of the polarization, the 
fully polarized Fock state in Eq. fl50| ) becomes Eq. (P5|) which is a generalized su(2) coherent 
stated 

1*0 =«^-> b», -*> = t x (.;:„) 1/2 e-i*n> (re, 

which becomes clear if one makes a correspondence between Eqs ( |7D|) and (fi5| ) as 

*- B ^ m , * + „^- m ), { ^ e *. or 
jo + w — ► , jo — n —> (M — m) , £ — > e 

where in Eq. (|70|) |jo^) = I jo — jo + w). Here, J± are the standard stt(2) raising and 
lowering operators of the su(2) angular momentum algebra defined by the generators Jj ,(i — 
0,...,3) 

J = (h\ Si + at, S 2 )/2 
J\ = (h\ Si — at ho) 12 

! , (72) 
J 2 = (a-i a-2 + S 2 Si)/2 

J 3 = (h\ a 2 — S 2 h\)/2i , 

where, considering that (a|, Si), (S 2 , a 2 ) represent two independent boson pairs, we have the 
standard su{2) algebra [Jj, Jj] = ie^k Jk , (i, j, k — 1, 2, 3). Here the central invariant of the 
algebra is J 2 = Jo(Jo + 1) = jo(jo + 1) where j = M/2. Here, M being an eigenvalue of J 
describes the total number of particles in the Fock state P5|). The uncertainty relations for 
j's are given by 

(A J) (AJj) > M ( J k ) , ^ fc = 1, 2, 3 . (73) 

The fully polarized Fock state is nothing but the generalized coherent state of the free field 
su(2) angular momentum algebra. Under certain conditions Eq ([70|) also coincides with the 
su(2) minimum uncertainty states0~il minimizing Eq.(|7^); and, this fact has been recently 
explored in the current literature in the context of su(2) interferometry.il~0 
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The idea of su(2) interferometry is to create interference between two arbitrary input fields 
by using passive and active lossless optical devices to measure the relative temporal phase 
between the fields. For this purpose the measured operators of the su(2) interferometry 
are defined as in Eqs ( |72"|) or they are related to Eqs ( [72]) by certain unitary transformations 
induced by the passive and active optical devices. These transformations on Eqs ( [72] ) [or 
the inverse transformations on the initial fields] can be engineered in such a way that the 
relative phase shift between the input fields can be measured by pure intensity measurements 
on the fields at the output ports of the interferometer.0 The principles of the quantum 
interferometry are thus based on a generalized operational scheme that is, in principle, very 
similar to the idea of the operational phase measurement presented in [2,3,6,9] as well as the 
present work. 

Now let us construct the uncertainty product for the general fully polarized Fock state 
|M)^ 0)7 and particularly focus our attention on the specific limit |M)^ at 7 = tt/4. The 
measured interferometric operators correspond, in the standard su(2) interferometry, to the 
expected values of the operators in Eqs([7^) or some linear superpositions of them in the 
initial state. For |M)^ 7 being the initial state, we have 

(Jo) = f 
(J 1 )=fcos2 7 

(^2) = y sin 27 cos 0o 

(J3) = t sin 27 sin 0o 

and 

(AJ0) 2 = 
(AJO 2 = M sin 2 2 7 

V ; 4 (75) 
(AJ 2 ) 2 = f (1- cos 2 0o sin 2 2 7 ) 

(AJ 3 ) 2 = f(l- sin 2 0o sin 2 2 7 ) . 
Using 7 = 7r/4 in Eqs. ( [75] ) and ( |73| ) we observe that \M)^ is an important state in the 



algebra defined by the operators in Eq. ([72]) . It is an su(2) coherent state [see Eq. ([70]) ] as 
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well as a minimum uncertainty (intelligent) state minimizing Eq. (|73|) for i ^ j = 1, 3; k = 2 
and z 7^ j = 1, 2; = 3. This can be explicitly seen by using Eqs. ([74]) and ([75]) in Eq. ([73]). 
Furthermore, when 7 = 7r/4, this result is independent from 0o; hence, a temporal shift 
in 0o does not change any of these properties. This implies that, if sw(2)-interferometric 
techniquesii are employed for \M)^ , the standard precision can be achieved in the mea- 
surement of the temporal phase@0 The precision in the phase measurement can be found 
from Eqs (f?4]) and (|75|) for the general case with |M)^ 0i7 as 



SMl) = ^ = 1 Vj^Wj. ( 76) 
\d(Ji)/d<l>o\ VM sin27cos0 o 



where 5(7) > 50o( 7r /4) = 1/vM which is the well-known minimum standard noise limit. 
Hence, theoretically, the maximum precision in the phase measurement can be achieved only 
at 7(7r/4) = 1 corresponding to 77 = 1. The basic idea being the extraction of the phase 
statistics from pure photon counting, the su(2) interferometry is in a close analogy to the 
operational measurement scheme. The operators Jj , (i — 0, 1, 2, 3) are the interferometric 
analogs of the operational ones Sj , (i — 0, 1, 2, 3) in Eq. (|7J), but there are also significant 
differences between them. Although the Jj's are the generators of the su{2) algebra, the 
Ej's all commute with each other and no useful uncertainty product similar to Eq. ( ]73| ) can 
be written for them. Now a legitimate question arises on how much the properties of the 
quantum state |M)^ 0i7 , as far as the Jj's are concerned, are preserved in the operational 
measurement scheme using the Sj operators. The main difference arising from the presence 
of the vacuum states in the Ej's as well as the operational scheme itself, it is nevertheless 
expected that for sufficiently strong fields the quantum operational measurement using the 
Sj operators should be consistent with Eq. (|73|). The deviations in the quantum operational 
measurement scheme from Eq. (]73| ) are expected when the initial field is sufficiently weak. 
Hence, by examining the uncertainty properties of |M)^ 0)7 , particularly near 7 = 7r/4, a per- 
fect ground to understand the influence of the operational scheme in the final measurement 
can be provided. 

We start the analysis of the uncertainty relations for |M)<^ 0)7 by examining the 7 dependence 
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of the measured D(9) and D((j>). The results are represented in Fig. 13 for M = 1,5 and 
for linear polarization in the range < 7 < tt/2. The figure indicates that, similarly to 
the results obtained for the fully polarized coherent state measurements, it is not possible 
to simultaneously minimize the fluctuations in the measurements of the 9 and related 
moments. The values of D{9) and D{(j>) in Fig. 13 corresponding to 7 = tt/A (i.e., rj = 1), 
7 ~ 0.37T (i.e., 77 = 0.5), 7 ~ 0.47T (i.e., 77 = 0.1), and 7 ~ 0A7n (i.e., 77 = 0.01) can also 
be read from the Figs. 11 and 12. Here 7 = 7r/4 has a special importance since this point 
corresponds to where |M)<^ 0j7 becomes a coherent as well as a minimum uncertainty state 
of the Jj's. As 7 is shifted away from 7r/4, particularly towards 7 = 0, 7r/2, one particular 
polarization mode starts dominating where the initial state gradually starts looking like a 
single mode Fock state. The single mode Fock limit is realized at 7 = 0,7r/2 for 9 related 
measurements [i.e., D(9) = and P{9) comprises a single delta function peak], and a 
maximally random fluctuations are observed in the <fi related measurements [i.e., D{4>) = 1 
and P(4>) = 1/271"]. Because of the fact that the interferometric operators do not commute 
with each other, it is not possible to find the interferometric analogs of the trigonometric 
operators Cg, Sg and C^, defined in Eqs. (|S|). This implies that, the interferometric analogs 
of D(9) and D{(j>) cannot be found by direct analogy and a comparison between the theory 
and the measurement is not possible for them. At this level, the only comparison with the 
theory can be made by examining the minimum uncertainty product for the Jj's and the 
ti's. 

Keeping 7 as the parameter, we now express Eq. ( ff3|) in |M)^, 0)7 in the form 

and find from Eqs. flT4|) and ( [75]) that 

W( 7 ) = ^(l-isin 2 2 7 ) (78) 

with the minimum uncertainty corresponding to U(tt/A) = 1/4. The operational analog of 
Eq ( |f7|) in terms of the Ej's can be found by direct inspection of Eqs (^-0) and (|72] ) as 
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^(7) = ^/^: (A^) 2 (79) 



( AE 2 ) 2 + (Af^ 2 
(S2) 2 + (S3) 2 

with all fluctuations in Eq. (|79| ) calculated within the operational scheme outlined in Sec. (II). 
In comparing Eq. (|7£|) with Eq. (|77|) the differences arising from the different normalizing 
factors of the transmission and the reflection coefficients in the d^s in Eq. (||) should also 
be accounted for. The result of the numerical calculations for U ov {^i) jlAipi) is presented in 
Fig. 14 as a function of M. We also observed that Eq. (|79| ) has no 0o dependence for all M 
and 7 [not shown in Fig. 14] which is consistent with the theoretical calculation in Eq. ([T^) . 
The Fig. 14 indicates that the photodetection in the operational scheme unavoidably creates 
an additional noise in the measurement such that the theoretical value of the uncertainty 
is not reached until the initial state has sufficiently large (i.e. 5 ~ M) number of photons. 
Here, as 7 varies, there is a compromise between the value of the measured uncertainty 
product for large M and the detector noise for small M. For instance, at 7 = it/ '4 (i.e., 
77 = 1), the measured uncertainty product approaches the theoretical minimum uncertainty 
for large M, although the detector noise is as large as 100% at the small M limit. On the 
other hand, as 7 deviates from tt/4, the measured uncertainty product is no more at the 
minimum for large M but the detector noise is smaller for small M. Hence, it appears that 
there is no global optimum value for 7. We thus conclude that 7 can be optimally fixed only 
depending on the individual observables chosen in the measurement (i.e., a result which we 
have also reached in the fully polarized coherent state example in Sec. II.A.l). 
In the theoretical interferometric calculations it is a common practice to neglect the influence 
of the photodetection. This is certainly a valid assumption if the initial field is sufficiently 
strong. On the other hand, we expect the additional higher bound on the uncertainty product 
to be a manifestation of any scheme based on photon counting in the weak-field limit arising 
from the quantum nature of the photodetection. Hence it is also natural to expect these 
effects to be observable in the su(2) interferometric measurements. This result indeed needs 
experimental verification, particularly, considering the advantage that certain schemes have 
been proposed for the generation of such quantum states as Eq. ( pf5| ) experimentally using 
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active non-linear processes.l 



III. DISCUSSION 

In this work, we focused our attention on the operational measurement scheme as applied 
to certain fully polarized quantum states particularly in the weak-field regime. We have 
shown that, similarly to NFM's operational phase measurement scheme it is possible to 
base the measurement of the state of polarization on pure photocount measurements, hence 
providing another example for operational approach. In particular, the measurement of 
the fluctuations of the temporal phase between the polarized field components is, not sur- 
prisingly, identical to the original work by NFM. The statistical behaviour of the Stokes 
parameters is investigated in terms of the trigonometric operators in Eqs. (|j) and the oper- 
ational counterparts of the quantum Stokes parameters of the polarized field are introduced 
in Eqs. (|?|). The application of the operational polarization measurement scheme is made 
to fully polarized quantum coherent as well as Fock states. With the purpose of extract- 
ing the detectors' influence on the measurement, the statistics of the measured fluctuations 
are examined and compared with the theoretical calculations. Our results confirm those of 
NFM's operational phase measurement scheme to conclude that the photon-counting pro- 
cess introduces additional noise in the final statistics particularly in the weak-field regimes. 
For sufficiently strong fields, the operational measurement scheme is consistent with those 
theoretical predictions in which the photo-detection effects are not included. 
The connection between the operational approach to the measurement of polarization and 
the su{2) interferometry is examined and the uncertainty principle is used as a means of 
analyzing the photodetection effects in the measurement where applications are made on 
the fully polarized Fock state. 

Operational approach to the measurement of phase has been investigated by D'Ariano and 
Parish in the context of quantum estimation theory^ which provides a unified formulation 
of the measurement process and the initial system under investigation. The quantum prob- 
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ability distribution of the iV-port homodyne detection in Eq. (|9]) is a specific example of the 
probability-operator-valued measure (POM) in the quantum estimation theory. An ideal 
quantum measurement is realized when the POM is based on an orthogonal and complete 
set of states comprising the eigenspace of the measured observable. Hence, depending on 
the nature of the measured observables of the initial state, finding an optimum detection 
scheme is the primary goal of a unified formulation of the measurement and the initial sys- 
tem. For the measured observable being phase related quantities, such an approach has not 
been idealized yet because of the fact that an orthogonal POM cannot be physically realized 
for the phase observable. With this in mind, one resorts to optimizing the phase measure- 
ment by a proper choice of the initial states as well as the parameters of the measuring 
system. At this point a connection is present between the primary goal of the quantum 
estimation theory and the attempts to surpass the standard noise limit by using interfero- 
metric transformations on the measuring system (or inverse transformations on the initial 
state). The optimal choice of these transformations using active as well as passive optical 
devices naturally depends on the initial state. Furthermore, it is also desirable (and under 
certain conditions it is strictly required) to have the transformed state conserve the basic 
features of the original state i.e., full polarization, coherent and minimum uncertainty state 
for Eq. (f|5|), the statistics of the fluctuations, etc. It is natural that for the initial state be- 
ing fully polarized, the full polarization itself is a strict condition that should be conserved 
by the transformations. On the other hand, since |M)^ is a minimum uncertainty state, 
the quantum statistics of the temporal phase and the fluctuations in corresponding de- 
pendent operators are coupled with those describing the fluctuations in the relative photon 
number n\ — fi2. Hence, the minimum uncertainty condition of |M)^ will most certainly be 
at stake after such transformations and this will change the quantum nature of the state [for 
instance, a rotation in the field space by 7 does not change the full polarization property 
but changes the minimum uncertainty relations, as it can be seen from Eq. ( |73|) and (|78|)1. 
For those states which are not the minimum uncertainty ones, this observation is still valid 
to a lesser extent. We nevertheless conclude that attempts to surpass the standard noise 
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limit for the fully polarized quantum states have to comply with a number of restrictions, 
which certainly renders it to be rather interesting problem. 
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Figure Captions: 



Fig.l The experimental setup to measure the classical and quantum Stokes parameters. 
Note that this setup is also able to measure all of the total of six Stokes parameters for a 
partially polarized light. 

Fig. 2 Ellipsometry for the fully polarized transverse electric field E in the tangent plane. 
Angular parameters are shown as defined in Eqs (f|) . 

Fig. 3 Measured probability distribution P(9) versus 9 for the fully polarized coherent state 
and for a) 77 = 1.0, b) 77 — 0.5 for the indicated average total photon numbers. 

Fig. 4 Measured probability distribution P(4>) versus <fi for the fully polarized coherent state 
and for a) 77 = 1.0, b) 77 = 0.5 for the indicated average total photon numbers. 

Fig. 5 Second-order fluctuations in the 9 related measurement for the fully polarized co- 
herent state and for the indicated values of 77. 

Fig. 6 Second-order fluctuations in the <fi related measurement for the fully polarized co- 
herent state and for the indicated values of 77. 

Fig. 7 Second-order fluctuations in the 9 and related measurements for the fully polarized 
coherent state as a function of the rotation parameter 7 for the indicated values of 77 and 
the relative temporal phase in the extreme weak-field limit (So) = 1. 

Fig. 8 Second-order fluctuations in the 9 and </> related measurements for the fully polarized 
coherent state as a function of the rotation parameter 7 for the indicated values of 77 and 
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for (S ) = 9. 



Fig. 9 Measured probability distribution P{6) versus 6 for the fully polarized Fock state 
for a) 7] — 1, b) 7] — 0.5 and for the indicated average total photon numbers. 

Fig. 10 Measured probability distribution P(<f>) versus for the fully polarized Fock state 
for a) 7] — 1, b) 7] — 0.5 and for the indicated for the average total photon numbers. 

Fig. 11 Second-order fluctuations in the 9 related measurements for the fully polarized Fock 
state as a function of the average total number of photons and for the indicated values of r\ 
(here we considered O = 0). 

Fig. 12 Second-order fluctuations in the related measurements for the fully polarized Fock 
state as a function of the total number of photons and for the indicated values of 77 (here we 
considered O = 0). 

Fig. 13 Second-order fluctuations in the 9 and related measurements for the fully po- 
larized Fock state as a function of the rotation parameter 7 for the indicated values of M 
and 0o- 

Fig. 14 Comparison between the theoretical uncertainty product and the measured one in 
the fully polarized Fock state for the indicated values of 77. 
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